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1. Introduction 

The recent discovery of a Higgs-like boson at the Large Hadron Collider (LHC) [Q] repre¬ 
sents an important milestone on the understanding of mass generation and electroweak symmetry 
breaking (EWSB). Nevertheless, an explanation for the observed distribution of fermion masses is 
lacking. Hierarchies in this pattern suggest that while third-generation masses may arise at the clas¬ 
sical level, quantum corrections are responsible for the origin of the remaining ones [^. Realistic 
models of radiative fermion mass and mixing generation have been built both in the frameworks 
of supersymmetric theories 0 and Grand Unified Theories (GUTs) [^. Fermion mass hierarchies 
may alternatively be explained by the presence of symmetries acting on flavor space and constrain¬ 
ing the structure of the Yukawa couplings. One can consider, for instance, the Froggatt-Nielsen 
scenario |^] where U(l) shaping symmetries are broken when some scalar fields S (flavons) ac¬ 
quire nonvanishing vacuum expectation values (VEVs), {S). The effective Yukawa structures thus 
generated are governed by powers of small parameters e {S)/M. 

Concerning the explanation of the neutrino flavour pattern, where the mass hierarchy is mild 
and the mixing large, the recently taken route is to implement discrete symmetries |0]. One possi¬ 
ble model-building approach in this context is to extend the field content by adding several scalar 
doublets, whose couplings to fermions depend on the broken or unbroken symmetries, and may or 
not be generated at the renormalizable level. The VEV and Yukawa coupling configurations will 
then determine the pattern of fermion masses. If a certain VEV configuration cannot accommodate 
the observed values of fermion masses at tree level, it is customary to dismiss it as not phenomeno¬ 
logically viable. This would occur if, for instance, a charged lepton e, is massless (or too light) at 
tree level in a multi-Higgs doublet model where only one of the doublets carries a VEV. In this case, 
a sizable coupling of these leptons with the nonzero-VEV Higgs can nevertheless be induced by 
quantum corrections, contributing to a mass after EWSB. In the present work, the possibility of 
radiatively generating charged-lepton masses due to such corrections is explored. Simple examples 
in two- and three-Higgs doublet models are considered, and the case of NHDMs extended with 
right-handed neutrinos is also briefly discussed. 

2. Radiative charged-lepton masses in NHDMs 

Consider the extension of the Standard Model (SM) with N Higgs doublets 
(2,1 /2). The Yukawa Eagrangian reads 



( 2 . 1 ) 


where qu, (-u denote quark and lepton doublets, dRi and CRj correspond to the right-handed 
singlets, and <i>a = i 02 ^*a = ^~'PaY■ The Yukawa matrices Y^ are general 3x3 complex 
matrices which can be diagonalized by the biunitary transformations 


VfYfUf = diag(yfi,yf2,3'f3), 

with y^j are real and positive. After EWSB, mass matrices are generated at tree level as 


( 2 . 2 ) 
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where Va = {<pa) is the VEV of <j)^. We restrict our attention to corrections to the charged-lepton 
Yukawa couplings, Y^. At the one-loop level, the renormalization group equations (RGBs) for 
these couplings read 

= = (2A) 

with /I and A the renormalization and reference energy scales, respectively. For the NHDM, the 
beta function at one loop is |||] 



+ ^b~r ^' 






(2.5) 


with ttg = —9g^/4 — 15g'^/4, where g and g' are the SU(2)vr and U(l)y gauge couplings and 




ab 


= 3Tr(Yf YU +Tr(Y^Y: 


d-sjd^\ 




(2.6) 


The contribution controlled by these coefficients stems from Higgs wave-function diagrams with 
ingoing <I>a, quarks and leptons in the loop, and outgoing <!>;,. 

In order to simplify our analysis, we henceforth consider the Higgs basis |^] (although the 
analysis is valid for a general vacuum configuration) where all VEVs are zero except the one of 
some 0Q (Va = V = 174 GeV), and Yukawa couplings are rotated accordingly. The charged-lepton 
mass matrix is therefore 

M, = vYf , VfYfuf = diag(y„y^,y,) (2.7) 


at tree level. Notice that we are not interested here in possible flavor-changing neufral currenf 
consfrainfs [10|, since fhese can be avoided considering fhe decoupling limif of NHDMs [|TT|]. 

We will focus in fhe one-loop correcfed couplings, denoted by Y^^^ The befa function for Y^ 
given in Eq. ( ^7^ confains terms depending on Y^ itself, contributing to existing tree-level masses, 
and terms proportional to Y^, namely the aforementioned oty^Y^ term. For b^a, such terms induce 
corrections to Y^ which are independent from its own structure. Under the reasonable assumption 
that the top quark couples with the same strength to all Higgs doublets, y, ~ 1, and considering, for 
simplicity, all matrices Y^ diagonal with elements given by y^i, one can estimate the magnitude of 
corrections to Y^ coming from Y^ to be typically of the order of 


5yai ~ 


^ybi 

1671 ^ 


log 



i = 1,2,3 , b ^ a, 


( 2 . 8 ) 


where A is a high scale at which the Yukawa couplings are initially given, and mn is the typical 
scale of the extra scalars in the theory. This illustrates the fact that Yukawa couplings with the 
zero-VEV scalars, while not contributing to tree-level masses, can induce important corrections to 
Y^, depending on A, mn and on the size of the Y^ [b ^ a). In order to end up with three massive 
charged leptons after taking quantum corrections into account, one requires ^(Y^^^) = 3, where r 
denotes matrix rank. It can be seen that, in general, the rank of Yukawa matrices in a NHDM may 
change due to these corrections, and lepton masses which were absent at tree level are induced. 
One starts by noticing that the beta function for Y^, at any loop order n, can be cast in the form 




k=\ 


(2.9) 
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where X^”] are matrices in flavor space. For instance, from Eq. ( |2.5[ ), one has at one loop 


Xl!] = {Sak + af) 1 + YfYf +15,, YfYf ^. 


( 2 . 10 ) 


In the leading-log approximation, the structure of the couplings Y^^”^ (corrected up to order n) will 
be given in the form of Eq. ( ^^ with different X matrices, In the general case where at least 
one {b 7 ^ a) is nonzero, maximum rank for these X' matrices is expected already at one-loop 


level. Since r(AB) < min{r(A),r(B)}, it is clear that the rank of each term in the sum of Eq. 
is at most equal to r(Y^), while since rank is subadditive, i.e. r(A-|-B) < r(A) -|-r(B). One thus 
has 


N 


N 


r{Y7 




. k=\ 


< minp,£ r(Yf 


( 2 . 11 ) 


k=\ 


Barring tailored cancellations among Yukawa structures, equality generally holds in Eq. ( |2.11] ). 
This result implies that one can in principle generate charged-lepton masses from low-rank tree- 
-level Yukawa couplings. Hence, the number of massive charged leptons will depend on these ranks 
and on the number of existing Higgs doublets. Eor instance, to end up with rim massive charged 
leptons from rank-1 (tree-level) mass and Yukawa matrices one would need^ N > n,„- To make this 
point clear, we consider in what follows some simple and extreme examples where radiative mass 
generation is key. It is not our goal to frame these examples in specific models but to provide a 
proof of concept for these claims. 


Considering first the 2HDM case Q12|], one can see that at most one charged-lepton mass can 
be radiatively generated from rank-1 Yukawa matrices in the Higgs basis. Instead, if 3 — n charged 
leptons are massive at tree level, one needs ^(Y^) > n to end up with all charged leptons massive. 

In the 3HDM, on the other hand, all three charged leptons could acquire a mass from rank-1 
Yukawa matrices in the Higgs basis. To illustrate this, consider a 3HDM with a vacuum config¬ 
uration of the type (vi,V 2 ,V 3 ) = (0,0,v) and with Yukawa couplings^ Y^ = diag(ei,0,0), Y^ = 
diag( 0 ,e 2 ) 0 ), Y 3 = diag( 0 , 0 ,e 3 ), taking all parameters real thus neglecting possible CP-violating 
effects in this sector [|^]. Within this setup, only the tau is massive at tree level and only Y 3 is rele¬ 
vant for charged-lepton masses. At one loop, in the leading-log approximation, the charged-lepton 
masses thus obtained read 


m, 


( 1 ) 


I6n^ 


tty 


e, log ( — ) , i= 1,2 . 
Mh/ 


( 2 . 12 ) 


Hereafter we consider A = Aqut ~ 10^® GeV and mh ~ 1 TeV. Keeping ay ~ '^(1)> the right 
masses are in this case obtained for (ei ,e 25 £ 3 ) — (1-6 x 10^^,3 .2 x 10^^,0.01), whose hierarchical 
structure could be explained by a Eroggatt-Nielsen mechanism. 

Even though so far only extreme examples of mass generation have been considered, less triv¬ 
ial Yukawa structures can also introduce corrections to lepton mixing. This effect can be illustrated 

'This condition is necessary but not sufficient due to the aforementioned possibility of cancellations, a trivial such 
exception being Yj. 

^Notice that in the following examples the zero entries in the Yukawa matrices should perhaps not be taken as 
strict zeros but interpreted as the limit of very suppressed entries. Our results are valid also for more complicated flavor 
structures which we do not consider here. 
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by considering the same 3HDM as before with a nondiagonal choice of rank-l matrices: 



( £1 0 0\ 
-£l£ 0 0 

.y5 = 

^0 £2£ 0^ 
0 £2 0 

II 

^0 0 £3£^ 
0 0 0 


-£l£ 0 Oy 


0 

0 

0 


y 0 0 £3 ! 


(2.13) 


As before, only one charged lepton is massive at tree level, while taking radiative effects into 
account corrects Y 3 to Y^^'^ = Y 3 — 5 Y 3 , with 

I '\\ ^2 ^ 


5Yi ~ 
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IhTT^ 


ay^Ei ape2e ^£3e(e|£^ + 6e-?' 
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-ay^£i£ ape2 
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-ay^£l£ 
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£|£3£^ 

3 3 

2^3 


log 


(-) 

\mHj 


(2.14) 


and tty given by Eq. ( |2.6| ). At this (one-loop) level, one obtains 

i;3l£j ,- / A \ /t^32 

IhTT^ 


„,i» ^yiWviog (Al. 4') ^yiTiiviog (Al 


(2.15) 


from the mass matrix ’ = VY 3 . Taking (£i,£2,£3) ~ (1.5 x 10^^,3.2 x 10^^,0.05) and £ <C 1, 
while keeping tty' ~ ^(l)^ allows one to reproduce the observed values for Moreover, the 

left-handed rotation which brings to the diagonal basis is no longer trivial: 

/-1 + £2 £ £ \ 

Vi.~ 


1- 


V 


0 


-£^ 1- 


(2.16) 
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Such a rotation would be invaluable to correct, for instance, a tribimaximal (TBM) mixing pat¬ 


tern [ |14{ ] arising at tree level, as its exact version is presently excluded due to a nonzero re¬ 
actor neutrino angle. Were this the case in the above example, one would have the correction 
Utbm 

fields to the diagonal basis. The new mixing angles would then be: 

1 


V^Utbm to the lepton mixing matrix due to the rotation of charged-lepton left-handed 


. 2 ^ 14 .2^ 

sin 012 ~ - + -£ , sm 623 


-(1 - £ ) , sin @13 ~ 2 £ 


(2.17) 


which, for £ ~ — 0 . 11 , lead to the right value for the reactor neutrino angle 613 , while keeping the 


remaining mixing angles within experimentally allowed ranges []15|]. This simple example shows 
how a scenario which would be excluded by tree-level considerations becomes phenomenologically 
viable when quantum corrections to the charged-lepton Yukawa couplings are included. We point 
out that we have started from a situation where the muon and the electron were massless and 
013 = 0, to end up with a case where and 613 are radiatively generated. Quark masses and 

mixing can also be generated through this same effect |[T^. 

Although we have presented examples with massless charged leptons at tree level—which 
would certainly call for a justification—we stress that, even if this is not the case, the RG corrections 
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to the Yukawas in NHDMs should always be kept in mind. To illustrate this, consider a 2HDM 
with real and diagonal Yukawa matrices Y^ = diag(yai ,ya 2 ,ya 3 ) in the Higgs basis where {(l)^) = 0. 
The one-loop corrected Yukawas will be 


( 1 ) 

y\i -yii 


{SltHg + ay’’] 

\6n^ 


yti^og 


(-) 

K^hJ 


(2.18) 


Taking the natural value aP' ~ ^(1), it is apparent from the above estimate that if y 2 i > \6n^y\il 
log(A/m//), the one-loop contributions to the Yukawa couplings are relevant. If, for instance, the 
electron couples with <I>i and <I >2 with yjg = m^/v ~ 2.9 x 10^^ and 4.6x10 /\og{K/mH), 
the correction to its coupling will be of the same order (or larger) than the tree-level value. Choosing 
as before A 10^® GeV and mn 1 TeV, this bound becomes y 2 e > 3.5 x 10 which may still 
be a relatively small parameter. For the third generation, y 2 T ^0.1 would lead to a correction of 
the order oiyu = m^/v. In fact, even if each charged lepton couples with the same strength with 
the two Higgs doublets, the Yukawa couplings would become = yn — Syu, where 


5yi, ^ {cCg -|- ccp cip) / A \ 
yii 16712 


(2.19) 


Taking Ug -f aP -I- ap ~ 4, A ~ 10^® GeV, and mh ~ 1 TeV, the induced dyn amounts to ~ 70% 
of the tree-level masses. 

It is worth pointing out that loop effects not controlled by the RGB running, though out of the 
scope of the present analysis, may need to be considered. These contributions have recently been 
shown to be at most 5% in several versions of the 2HDM []T7[]. 


3. Radiative charged-lepton masses in NHDMs with right-handed neutrinos 


Consider now the possibility of extending an NHDM by adding three right-handed neutrinos 
VRj to the field content, with masses Mj » v. The extended Yukawa Lagrangian reads 


^' = ^- [{YDijhi^aVRj + H.C.] , 


(3.1) 


where ^ is that of Eq. ( |2.1[ ), while Dirac neutrino Yukawa matrices are denoted by Y^. 

It is known that in the SM the presence of couplings Y'' cannot generate charged-lepton 
masses radiatively. This is due to the fact that new terms in the beta function of Y^ are of the 
form Y''Y''^Y^, implying that Yukawa eigenvalues are always proportional to themselves |18]. In 
contrast, different Y^ in the NHDM will contribute to the renormalization of the various Y^, both 
through wave-function and vertex corrections above the right-handed neutrino mass scale. It can 
be shown that the modified one-loop beta function, valid from A > M, to M,, reads 


K'" = Pi" + « + f Y^Y'tYj - 2YJYyY 


b 1 


(3.2) 


with af = Tr(Yf Y^), in agreement with Ref. Ql^. The presence of right-handed neutrinos adds 
two new contributions to the beta function which do not depend on Y^ [cf. Eqs. ( |2.5[ ) and (3.2)]. 
Notice, however, that these extra terms are only active for /r € [Mi, A > Mi]. 
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It can be seen that the second term of Eq. (3.2) alone can generate masses for both the electron 
and the muon in an NHDM. Extending the previously presented 3HDM example by the addition 
of right-handed neutrinos, and considering, for simplicity, diagonal Yukawa matrices in the Higgs 
basis, one obtains 


m, 


( 1 ) 


^ 21 


eilog 


A 

Mi 


/ = 1 , 2 , 


(3.3) 


where e, are the couplings of the electron (/ = 1) and the muon (/ = 2) to the zero-VEV Higgs 
doublet <I>i. Taking YJ ~ ^(1), which implies ~ seesaw (light) neutrino masses my ~ 
0.05 eV require a heavy scale of M, ~ 10^"^ GeV. In this case, one could obtain the right electron and 
muon masses for (£i,£ 2 ) — (10^^,0.02). If, on the other hand, one assumes that charged leptons 
couple with the same strength to the doublets <I>i and <I> 2 , the corrections to the masses coming 
from the Yj’s do not exceed 10% of the tree-level values. 


4. Conclusions 

In this work it has been shown that charged-lepton masses can be radiatively generated in a 
natural way in NHDMs due to the presence of Yukawa couplings of leptons and quarks with the 
extra (non-SM) Higgs doublets. It is important to mention that, besides affecting charged-lepton 
masses, the RGE running may have dramatic effects on the left-handed rotation which brings the 
charged-lepton mass matrix to its diagonal form. This clearly changes predictions for lepton mixing 
parameters, which have to be confronted with neutrino oscillation data. These corrections should 
then be taken into account in phenomenological studies of fermion mass and mixing models with 
more than one Higgs doublet. The presence of right-handed neutrinos provides additional sources 
for this effect. 
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